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Abstract 

The theory of gamma-ray production in solar flares is treated in detail. 
Both lines and continuum are produced. The strongest line predicted at 
2.225 MeV with a width of less than 100 eV and detected at 2.24 + 0.02 MeV, 
is due to neutron capture by protons in the photosphere. Its intensity is 
dependent on the photospheric 3 He abundance. The neutrons are produced in 
nuclear reactions of flare accelerated particles which also produce positrons 
and prompt nuclear deexcitation lines. The strongest prompt lines are at 
4.43 MeV from 12 C and at ~6.2 from 16 0 and 16 N. These lines result from both 
direct excitation and spallation. The widths of individual prompt lines are 
determined by nuclear kinematics. The width of the 4.43 MeV line is ~100 keV 
and that of the 6.3 MeV feature is ~300 keV. Both these lines have been 
observed from a solar flare. Other potentially observable lines are predicted 
at 0.845 and 1.24 MeV from 06 Fe, at 1.63 MeV principally from 14 N and 20 Ne, 
at 1.78 MeV from ss Si, at ~5.3 MeV from 1E 0 and 15 N, and at 7.12 MeV from 16 0. 
The widths of the iron lines are only a few keV, while those of the other lines 
are about 100 keV. The only other observed line is at 0.511 MeV from positron 
annihilation. The width of this line is determined by the temperature, and its 
temporal variation depends on the density of the ambient medium in the anni- 
hilation region. Positrons can also annihilate from the 3 S state of positronium 
to produce a 3-photon continuum below 0.511 MeV. In addition, the lines of 
7 Li and 7 Be at 0.478 keV and 0.431 keV, which have kinematical widths of 
~30 keV, blend into a strong feature just below the 0.511 MeV line. 


From the comparison of the observed and calculated intensities of the 
line at 4.4 MeV to that of the 2.2 MeV line it is possible to obtain infor- 
mation on the spectrum of accelerated nuclei in flares. Moreover, from the 
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absolute intensities of these lines the total number of accelerated nuclei 
at the Sun and their heating of the flare region can be estimated. We find 
that about 10 33 protons of energies greater than 30 MeV were produced in the 
1972, August 4 flare. 

The gamma-ray continuum, produced by electron bremsstrahlung, allows the 
determination of the spectrum and number of accelerated electrons in the MeV 
region. From the comparison of the line and continuum intensities we find a 
proton-to-electron ratio of about 10 to 10 s at the same energy for the 
1972, August 4 flare. For the same flare the protons above 2.5 MeV which are 
responsible for the gamma-ray emission produce a few percent of the heat 
generated by the electrons which make the hard x-rays above 20 keV. 
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INTRODUCTION 

Measurements of accelerated charged particles near the Earth clearly 
indicate that such particles are produced in great profusion in solar 
flares. These particles consist of both electrons and nuclei. But until 
the advent of solar gamma-ray astronomy, observations in the radio and 
x-ray bands had revealed only the existence of the electronic component 
in the flare region itself. In the hopes of finding some measure of the 
properties of accelerated protons and heavier nuclei in flares, a variety 
of theoretical studies of the possible nuclear reactions of such particles 
in the flare region have been made by Biermann et al. (1951), Morrison 
(1958), Lapointe (I960), Hess (1962), Chupp (1964), Lingenfelter et al. 
(1965a, b), Dolan and Fazio (1965), Bland (1966), Lingenfelter and Ramaty 
(1967), Ito et al. (1968), Kuzhevskii (1968), Ito and Okazoe (1969), 
Lingenfelter (1969) and Cheng (1972). 

In a previous paper (Lingenfelter and Ramaty 1967) we treated in 
considerable detail the nuclear reactions produced by accelerated charged 
particles in solar flares and we showed that in large flares these 
reactions produce detectable lines in the gamma -ray region* We found 
that the strongest lines should be at 0.5, 2*2, 4*4, and 6*1 MeV resulting 
from position annihilation, neutron capture on hydrogen, and deexcitation 
of excited states in 12 C and 16 0, respectively. 

The recent observations by Chupp et al. (1973) of the first gamma- 
ray lines from solar flares confirm these predictions. During the flash 
phase of the 1972, August 4 flare all of these lines were observed with 
relative intensities essentially consistent with our calculations. 
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Since these observations became available, several additional studies 
on gamma-ray line production in solar flares have been undertaken: 

Ramaty and Lingenfelter (1973a) have investigated the consistency of 
the observations with the theory of nuclear reactions in flares; these 
authors also considered the effects of positronium formation and neutron 
propagation in the solar atmosphere (Ramaty and Lingenfelter 1973b) ; 
Reppin et al. (1973) treated the time dependence of the 2,2 MeV line: 

Wang and Ramaty (1974) have done a detailed calculation on neutron propa- 
gation and 2o2 MeV line formation and they pointed out the importance of 
photospheric 3 He a s a nonradiative sink for the neutrons; Kozlovsky and 
Ramaty (1974a) have pointed out that otot reactions produce the 7 Li and 
7 Be lines at 478 keV and 431 keV which could be observable from flares; 
and Suri et al. (1975) have measured the continuum in the gamma -ray 
region for the 1972, August 4 flare and discussed its significance in 
terms of two stages of particle acceleration* 

In the present paper we wish to summarize the above material and 
to present updated calculations on the production of gamma rays in 
solar flares. In Section 2 we define the interaction models that we 
use in our calculations; in Section 3 we consider neutron production and 
2o2 MeV line formation; in Section 4 we consider the production of prompt 
gamma-ray lines; in Section 5 we treat the continuum emission; in Section 
6 we compare the results of Section 3, 4 and 5 with data for the 1972, 
August 4 flare and we deduce the number and spectrum of accelerated 
particles at the Sun and the energy that they deposit in the solar 
atmosphere; in Section 7 we treat problems concerning the formation of 



the 0,51 MeV line; in Section 8 we provide estimates for the expected 
high-energy gamma -ray and neutron fluxes at Earth; and we summarize our 
results in Section 9* Some of the material contained in the present 
paper has been recently given in two symposium proceedings (Ramaty and 
Lingenfelter 1975a ,b) c 
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INTERACTION MODELS 

We consider an interaction volume in which a population of acceler- 
ated particles interacts with the ambient medium assumed to be 
uniform throughout this volume. For the composition of this medium we 
use the abundances of Cameron (1973) „ The instantaneous production, rate 
of secondary particles is given by 

q(E s ,t) = nj dE N(E t t) eg ff (E) f(E,E fl ) , (1) 

where n is the number density of the ambient medium, N(E,t) is the 
instantaneous number of accelerated particles per unit energy per nucleon 
at time t; E and E s are the energies per nucleon of the primary and 
secondary particles; eg is the velocity of the primary particles*, <j(E) is 
the cross section for the production of secondaries as a function of E? 
and f(E,E s )dE s is the probability that a secondary particle produced by 
a primary particle of energy per nucleon E will have energy per nucleon 
in dE g around E s . 

If we assume that the primary particles are produced in the inter- 
action volume at a rate N(E»t), and if they can only lose energy or 
escape from this volume, then the instantaneous number N(E t t) is related 
to the source function fr(E,t) by the continuity equation 

+ r| [ E (E , t ) N(E,t)] + I N (E , t ) = N (E,t), (2) 

where £ is the rate of energy loss and is the mean escape time. Both 
these functions can depend on energy and time, but we shall assume that 
they are only energy dependent « The solution of equation (2) can then 
be written as 

E 

= E" 1 f dE' N [ E 1 , t ' ] exp [ - T 1^ dE” T, 

t E T E(E”) J 


N(E,t) 


(3) 



where 


t 


dE" 

fi(E") 


(4) 


t' 


{ 


(e.g* Syrova tkii 1959, Lingenfelter and Ramaty 1967) 

Equation (3) can be reduced to simpler forms in two limiting cases * 
In the thin-target model, the particles escape rapidly from the inter- 
action volume * For this model to be valid, the escape time has to be 
sufficiently short so that N does not vary appreciably over a time 
interval of order t a and an energy interval from E* to E, where E'-E 
is the energy lost during t 1 . Also TE must remain essentially constant 
from E 1 to E* Equation (3) can then be written as 

N(E,t) = t x (E) N(E,t). (3) 


In the thick-target model the particles lose energy and stop in 
the interaction volume, i .e . t » t - t 1 at all energies of interest. 
Equation (3) then simplifies to 

N(E,t) = E-i I dE' N (E 1 , t ' ) , (6) 

where t T is given by equation (4)o 

We can now substitute equations (5) or (6) into equation (1). For 
the thin- target model we obtain 

q(E s , t) = n £ <JE cj3a (E) f (E,E g ) ^ (E) N (E,t), (7) 

while for the thick-target model we get 

GO . OD 

q(E g ,t) = n J dE cfcr (E) f (E,E g ) E" 1 J dE * N (E\t'). (8) 

E 

Since equation (7) is of the same form as equation (1) which is the more 
general of the two, we shall use the latter in our subsequent calculations 
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for thin-targets « If the conditions for a thin target are fulfilled, we 
can relate, by equation (5), the instantaneous spectrum N(E) to the 
source spectrum ft(E) . But even if these conditions are not satisfied, 
we can still use the thin-target calculations provided that N is inter- 
preted only as the instantaneous number of accelerated particles in the 
interaction region. 

Equation (8) can be further simplified by considering the time- 
integrated production rate 

Q(E S ) = r” q(E ,t) = n J dEcp a (E)f (E,E S ) IT 1 ]* dE'N(E'), (9) 

-*> * E 

where 

CO 

N(E') = J dt' M (E'.t 1 ). (10) 

-oo 

Because the inner integral in equation (9) is just the integral spectrum 
of tf(E), we have 

Q (E s ) = A [“ dE dx <j(E)f(E,E s ) fT(>E) . (11) 

* dE 

where d£ _ EA/(ncp) , and A is Avogadro's number. If we are not inter- 
ested in the energy of the secondary particle we can integrate both 
equations (1) and (11) over E s , This yields 

CO 

q = n | dE cpcj(E) N(E) (12) 


and 


Q = A J* dE dx a (E) N (>E) 
• dE 


(13) 


We shall use equations (12) and (13) in our subsequent calculations 
of nuclear reactions where we do not need the energy of the secondary 
particles. In the thin-target formula (12), N(E) is the instantaneous 
number of accelerated particles in the interaction region per unit 
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energy per nucleon at E. In the thick-target formula (13), N (>E) is 
the total number of accelerated particles of energies per nucleon 
greater than E that stop in the interaction region. For both N(E) and 
N (>E) we use power-law and exponential spectra* These are 

N ± (E) or N^(E) = (k ± or k ± ) E _S (14) 

and 

N ± (>E) or N ± (>E) = (k^ or k^) exp( -P/Po) . (15) 

Here i denotes the chemical or isotopic species of the particles; P is 
particle rigidity; the k f s are constants determined by normalizing the 
particle numbers to 1 proton of energy greater than 30 MeV, using the 
composition of the ambient solar atmosphere; and s and Po are, respectively, 
the spectral index and characteristic rigidity assumed to be the same for 
all charged particle components. 

Before proceeding with the calculations, we wish to comment on the 
assumed abundances for both the accelerated particles and the ambient 
medium. As we shall see, the prompt deexcitation lines are produced mainly 
by reactions induced by protons. Hence their intensities relative to each 
other depend on the relative abundances of nuclei with A^12 (mainly 12 C and 
16 0), and these are reasonably well known in the solar atmosphere. Neutron 
production, however, depends strongly on the helium abundance. From the 
data of Cameron (1973) we have that He/H « 0.069 and He/(C+N+0) =60. If 
the abundance of He relative to both H and CN0 is larger, then more neutrons, 
and 2.2 MeV photons, are produced relative to the prompt lines than indicated 
by our calculations. In particular if He/H = 0.1, the intensity of the 
2.2 MeV line relative to that of the 4.4 MeV line will be increased by 


about 40%. 
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Neutron and 2.23 MeV Gamma Ray Production 

Neutron production by accelerated charged particles was treated 
by Lingenfelter et al. (1965a) and Lingenfelter and Ramaty (1967) • In 
the present paper we have updated the cross sections used by these 
authors and we have added data at low energies including ty~particle 
induced reactions. 

The most important neutron producing reactions and their 
threshold energies are listed in Table 1. The cross sections for reaction 
(1) are the same as given by Lingenfelter and Ramaty (1967). Below the 
multiple TT production threshold these cross sections are obtained by 
subtracting from the total proton-proton inelastic cross section 
(Hess 1958) the cross section for the reaction 1 H(p, TT + ) s H (see Ramaty 
and Lingenfelter 1969). At higher energies, the neutron production 


cross section is assumed to be 


inelas 


, . v J-liCXd£> 

a pp “ (2 - kp) (jpp , (lb) 

ine las 

where c is the total inelastic cross section, approximately 

equal to a constant at 27 mb at these energies, and k^ is the proton 
multiplicity having approximately a constant value of 1.25 (see Lingen- 
felter and Ramaty 1967). 

The cross sections for the reactions (2) are from Meyer (1972). 

Above a few hundred MeV, these reactions include also the formation of 
pionSo The cross sections for reactions (3), (4), (5) and (6) above about 
15 MeV are from Lingenfelter and Ramaty (1967). At lower energies most 
of the neutrons from proton induced processes are from the reactions 


13 C (p,n) 13 N, 18 0(p,n) 18 F, 14 N(p,n) 14 0, and 5S Fe(p,n) 56 Co which have very 
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low threshold energies. The cross sections for these reactions are 
from Gibbons and "Macklin (1959), Bair (1973), Audouze et al. (1967), 
and Tanaka and Furukawa (1959). The cross sections for the fusion 
reaction (8) were summarized by Kozlovsky and Ramaty (1974b). 

The cross section for (a,n) reactions on l2 C, l3 C, 18 0, 33 Ne, 35 Mg, 
3S Mg and 39 Si are known experimentally at low energies. The cross 
section for 12c (a,n) l5 0 is from Nelson et al. (1963); the cross section 
for l3 C (a,n) 16 0 was taken from Walton et al. (1957), Sekharan et al. 
(1967), and Bair and Haas (1973); the cross section for 18 0(<y,n) 31 Ne is 
from Bair and Willard (1962) and Bair and Haas (1973); the cross section 
for 23 Ne(a,n) 3e Mg is from As he ry (1969) and from Haas and Bair (1973); 
the cross section for 2B Mg(a,n) 29 Si is from feair and Willard (1962) and 
we assumed the same cross section for 36 Mg (a,n) 28 Si; the cross section 
for 29 Si(a,n) 33 S is from Gibbons and Macklin (1959). There are no data 
for the reactions 14 N(a,n) l7 F, 16 0(« ,n) 19 Ne, and 30 Ne ( a ,n) 23 Mg. By 
assuming that these reactions have the same cross sections as that of 
the reaction 13 C(c*,n) ls O (Nelson et al. 1963), we find that (a,n) 
reactions on 14 N and ie 0 do not contribute significantly to neutron 
production in any energy region. The reaction ls O(a,pn) ls F for which the 
cross section was measured by Furukawa and Tanaka (1961) does contri- 
bute significantly to neutron production by incident ^-particles of 
about 10 MeV /nucleon. 

The cross section for the reaction Fe(^,n) was estimated theo- 
retically. Above the threshold for this reaction the emission of 
neutrons becomes the dominant reaction and therefore o(a» n ) ~ a c (e?) 
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where o c (ct) is the cross section for forming the compound nucleus with 
an alpha particle. For Fe the value of o c (a) in the few MeV/nucleon 
range can be obtained by using the continuum theory (Blatt and Weisskopf 
1952). Our numerical values are from Blatt and Weisskopf T s Table 4,1 
(p.353), averaged for = 1.3 x 1CT 13 and 1.5 x l(T l3 cm. 

The resultant cross sections are given in Figure 1 as functions of 
energy per nucleon E. Here pp, p&, pCNO and $CNO indicate neutron 

production in proton-hydrogen, proton-helium, alpha particle-helium, 
proton-heavy nuclei, and a part icle -heavy nuclei reactions, respectively. 

The pCNO and aCNO cross section are defined as 

6^cno = £ G''. _Ti (17) 

i 



where n. is the abundance of isotope i; o’ - and a are the cross sections 

pi & i 

for neutron production by protons and alpha particles on isotope i? and 
the summations are over all isotopes equal to or heavier than 
1S C. The cross sections for the inverse reactions 0P, CNOp and CNO& are 
the same as those for the direct reactions at the same energy per nucleon. 

The instantaneous neutron production rates in the thin-target model 
for power-law and exponential spectra are shown in Figures 2 and 3, 
respectively , The various production modes are: pp (proton-hydrogen), 
pot (proton-helium), (# particle-hydrogen), got (a part icle -he lium) , 
pCNO (proton-heavy nuclei), CNOp (heavy nuclei-hydrogen) , <*CN0 (a 
particle-heavy nuclei), and CNOa (heavy nuclei -he lium) . As can be seen 
for flatter spectra (smaller values of s or larger values of P Q ) the 
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neutrons are produced mainly in p$, pp and reactions. For power-law 
spectra (Figure 2) neutron production at large values of s is mainly 
due to ^CNO and CNOq> reactions; the contribution of pCNO and CNOp 
reactions is small at all values of s ; and reactions make a major 
contribution around s=4„ For exponential spectra, (Figure 3) almost 
all of the neutrons are produced in p# reactions at values of Po>25 MV. 

In this case the relative contributions of reactions induced by ^-particles 
and heavy nuclei are lower than for power-law spectra. Because particles 
with Z 2 have larger rigidities than protons, it follows from equation 
(15) that their fluxes relative to the proton flux at the same energy 
per nucleon is lower than in the power-law case given by equation (14). 
Since the nuclear cross sections are the same for the direct and inverse 
reactions at the same energy per nucleon, for particle spectra which are 
exponential in rigidity more neutrons are produced by a proton induced 
reaction than by the corresponding inverse process * 

The total neutron production rates in the thin-target model, q, 
and the total neutron yields in the thick-target model, Q, are shown in 
Figure 4 for power-law and exponential spectra . In the thin-target 
model, for flat primary spectra the neutron production rates are about 
the same for the exponential and power-law cases 0 This result is due 
simply to the fact that for such spectra most of the neutrons are pro- 
duced in p£> reactions with effective threshold around 30 MeV/nucleon and 
all of the assumed spectra are normalized of 1 proton above this energy. 

But because a power- law spectrum contains a much larger number of low- 
energy particles than an exponential spectrum with this same normalization, 
the steep power— law spectrum can yield orders of magnitude more neutrons 
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from ^CNO and CNOa reactions which have thresholds more than an order of 
magnitude below 30 MeV/nucleon. 

Similar effects are evident also in the thick-target model except 
that here the relative contribution of the low-energy particles in 
general is diminished because of their shorter range. 

Having considered the production of neutrons, let us now discuss 
their propagation and the ensuing gamma-ray line production. Wang and 
Ramaty (1974) considered in detail the effects of neutron propagation in 
the solar atmosphere on the production of gamma rays by the reaction 

n + p -» d + 7 . (19) 

In their treatment a distribution of neutrons was released in the 
chromosphere or corona, and the path of each neutron after its release 
was followed by a computer Monte-Carlo simulation. If the neutrons are 
released above the photosphere, any initially upward moving neutron 
escapes from the Sun. Some of the downward moving neutrons can also 
escape after being baekscattered elastically by ambient protons, but most 
of these neutrons either are captured or decay at the Sun. Because the 
probability for elastic scattering is much larger than the capture 
probability, the majority of the neutrons are thermalized before they 
get captured. Since the thermal speed in the photosphere (where most 
of the captures take place) is much smaller than the speed of light, 
the gamma-rays from reaction (19) are essentially at at 2.225 MeV and 
the Doppler-broadened width of this line is negligible (^ 100 eV) . 

The bulk of neutrons at the Sun are captured either on H or on 3 He. 
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Whereas capture on H yields a 2.2 MeV photon, capture on 3 He proceeds 
via the radiationless transition 

n + 3 He -*• 3 H + p, (20) 

and hence produces no photons. The cross sections for reactions (19) 
and (20) are 2.2 x 10 -30 e' 1 c# and 3.7 x 10" 2S p~ 1 cm 3 , respectively, 
where Cp is the velocity of the neutron (for details see Wang and 
Ramaty 1974). Thus^ if the 3 He/H ratio in the photosphere is ~ 5 x 10 -5 
comparable to that observed in the solar wind, nearly equal numbers of 
neutrons are captured on 3 He as on H. 

The results of the Monte-Carlo calculations of Wang and Ramaty 
(1974) are presented in Figures 5 and 6 for two assumptions on the 
photospheric ®He abundance: 3 He/H = 0 and 3 He/H = 5 x 10" 6 . In these 

calculations an isotropic distribution of monoenergetic neutrons of 
energy E n is released above the photosphere. The solid lines are the 
probabilities for the various indicated processes. As can be seen, the 
capture and loss probabilities increase with increasing energy, because 
higher energy neutrons penetrate deeper into the photosphere. This 
reduces their escape probability and leads to a shorter capture time, 
thereby reducing the decay probability. When 3 He/H = 5 x 10" 5 , the 
probability for loss on 3 He almost equals the capture probability on 
protons. The escape probability is greater than 0.5, because all initially 
upward moving neutrons escape from the Sun. Note- that the sum of all 
probabilities equals 1. 

The dashed lines in Figures 5 and 6 are photon yields per neutron, 
f(0>E n ), as a function of E r for various gamma-ray emission angles, 0, between the 
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earth-sun line and the vertical to the solar surface . 

At low neutron energies and e near zero, f is close to the capture 
probability on protons* This means that gamma rays from low-energy 
neutrons observed close to the vertical escape essentially unattenuated 
from the Sun* At higher energies and at larger angles, however, there 
is significant attenuation of the gamma rays due to Compton scattering 
in the photosphere. Even though f does depend on 0, for flares sufficiently 
close to longitude and latitude zero on the Sun we can neglect its angular 
dependence. This approximation is quite valid for the flare of 1972, 

August 4 since its longitude and latitude were EOS and N14. 

Using Figures 5 and 6, Wang (1975) has evaluated the average photon yield 
per neutron, f, by integrating f(E n ) over the neutron production energy 
spectra of the various interaction models discussed above. Such neutron 
spectra were evaluated by Lingenfelter and Ramaty (1967). Recently, 

Wang (1975) has revised these calculations, and some of his results are 
shown in Figure 7 and 8. The values of f obtained by using these spectra 
are shown in Figure 9. The upper and lower 4 curves are for 3 He/H = 0 
and 3 He/H = 5 x 1CT 5 , respectively. In terms of these average photon 
yields per neutron, the average 2.2 MeV photon flux at Earth in the 
thin- target model, 0, and the total 2.2 MeV photon flux 0 in the thick- 
target model are given by 


^2.2 

= 

qf/(4nR 2 ) 

(21) 

^2.2 

= 

Q f/C^R 2 ) 

(22) 


where R = 1 A.U. 
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We finally note that equation (21) is valid for the average neutron 
flux only, because the instantaneous 2„2 MeV flux lags behind the 
instantaneous neutron production rate. This lag is almost entirely due 
to the finite neutron capture time in the photosphere. Wang and 
Ramaty (1974) have investigated this effect, and some of their results 
are given in Table 2. Here <n> is the most probable density in the 
photosphere where the captures take place, t is the mean capture time, 
and is the neutron decay mean life. In terms of the parameter 
X = T c _1 + Td” 1 * the time profile of the 2.2 MeV photon flux from a 
monoenergetic burst of neutrons released at t Q can be approximated 
(Wang and Ramaty 1974) by 


0 (2.2 MeV) oc exp[-x(t-t Q )]. 


(23) 
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4. Prompt Gamma Ray 1 Lines 

The various prompt gamma-ray lines that can be produced in solar 
flares are listed! in Table 3. These lines result from the deexcitation 
of nuclear levels in isotopes which either are abundant in the solar 
atmosphere or are easily produced in nuclear reactions of energetic 
protons and alpha particles with the abundant constituents. Since 
1 H, 3 He and 4 He have no excited states which decay mainly by y-emission 
(Meyerhof and Tombrello 1968, Kozlovsky and Ramaty 1974c), the first 
group of nuclei consists of the isotopes 12 C, 14 N, 1 S 0, 20 Ne, 34 Mg, 28 Si, 56 Fe 
the second group is comprised of isotopes such as 6 Li, 7 Li, 7 Be, 10 B, 

1 :l B, 1 x C, l 3 C, l 5 N. In order to determine the strongest lines which 
might be expected, we considered laboratory experiments in which line 
emission was directly observed from proton and 3 -particle bombardment of 
targets which are abundant in the solar atmosphere. Such experiments 
were done by Clegg et al. (1.961), Foley et al. (1962), and Zobel et al. 

(1965, 1968) for 13 C, l 4 N, 16 0, 24 Mg and sS Fe targets, and incident 
protons and 3 -particles in the energy range from about 10 MeV to 150 MeV 
which essentially includes all solar flare particle energies of interest. 

All the lines in Table 3, except the lines at 0.431 and 0.478 MeV, have 
been observed in the laboratory. These latter two lines have been inferred, 
theoretically from the existence of bound states in 7 Li and 7 Be which 
decay by y-emission and which can be produced in 33 reactions (Kozlovsky 
and Ramaty 1974a) „ 

The cross sections for the reactions l 3 C (p,p 1 ) 12 C * 4 043 , 12 C(a, 3 , ) 12 C * 4 * 43 , 


16 0(p,p , ) 16 0 * 6 ' 14 , and 16 0 ( 3 , 3 ' ) 16 0 * 6 014 are given in Figure 10. The 
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cross sections for the proton induced reactions were summarized by 
Lingenfelter and Ramaty (1967). The cross sections for the reaction 
12 C(a,«') l2 C* 4 ' 4 3 were taken from the measurements of Correlli et al. 
(1959), Yavin and Farwell (1959), Mitchell et al. (1964), Morgan and 
Hobbie (1970), Fawzi (1972) and Smith et al. (1973). The cross sections 
for the reaction 16 0(a ,ot ’ ) 1S 0* 6 ° 14 were taken from Correlli et al. 

(1959), Yavin and Farwell (1959), Harvey et al. (1964), Blatchley 
and Bent (1965), and Bergman and Hobbie (1971). 

Excited levels in nuclei can be populated not only by direct 
excitation reactions such as those discussed above, but also by spallation 
reactions in which protons and o'-particles break up heavier nuclei into 
lighter fragments that emerge from the reaction in excited states. 

In particular, the level ^C* 4 * 43 can be populated by the reaction 
16 0(p,p<y) 13 C* 4 * 43 . The cross section for this reaction, measured 
by Zobel et al. (1968), is given by the dashed line in Figure 10. As 
can be seen, above about 30 MeV this cross section is greater than the 
direct excitation cross section. Because in the solar atmosphere 16 0 is 
about twice as abundant as 12 C, spallation of 16 0 makes a very import- 
ant contribution to 4.43 MeV photon production. 

By using these cross sections and the interaction models discussed 
above, we have calculated the production rates of 13 C* and 16 0*. We must 
distinguish, however, between reactions induced by accelerated protons 
or ^-particles, and reactions induced by accelerated heavy nuclei. For 
the former, the Doppler widths of the lines are small (~ 100 keV), but 
for the latter the lines are significantly broadened by the motion of 
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the excited fast nucleus which has lost little kinetic energy in the 
interaction*. Because these lines are so broad that they cannot be 
resolved from the background with presently available instrumentation, 
in our treatment we consider the line intensities from proton and 
^-particle induced reactions only 0 

The production rates of 12 C* 4 ° 43 and l6 0'* 6 * 14 in the thin-target 
model are shown in Figures 11 and 12 for power law and exponential 
spectra, respectively* We see that for flat spectra, these lines are 
produced mainly by proton induced reactions, whereas for steep spectra 
(large values of s), the contributions of the ^-particles becomes 
important* For exponential spectra, proton-induced reactions are the 
principle source of the excited states at all values of P G ^ 2b MV. 

The l2 C* 4-43 production due to spallation reactions are also shown in 
Figures 11 and 12. As seen, these reactions make an important contri- 
bution to the total production for essentially all spectra except steep 
power laws for which the low-threshold direct excitation reactions 
are the main source of the excited states* 

The ratio of the total production of the level ^C* 4 ’ 43 to the 
neutron production is shown in Figure 13 for both the thin and thick- 
target models as a function of s , for power-law spectra, and P , for 
exponential spectra* As can be seen, for exponential spectra in both 
the thin and thick-target models 13 C* 4o43 /n decreases with increasing P Q * 
This results from the increase of the neutron production cross section 
with increasing energy as opposed to the decrease of the excitation 
cross sections (compare Figures 1 and 10) * The same behavior can be 
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seen for power- law spectra for s smaller than about 4.5 in the thin- 
target model, and s less than about 6 in the thick- target model. For 
larger values of s, 12 C* 4 * 43 /n decreases with increasing s because these 
neutrons are produced mainly in cy CNO and CNOcx reactions which have lower 
thresholds than those for prompt gamma-ray production. 

The 4.4 and 6.2 MeV lines are the most intense prompt lines and the 
only prompt lines so far observed for a solar flare. In Table 3 we list 
the approximate expected intensities relative to the line at 4.43 MeV 
of the various other gamma-ray lines which might be observed during flares. 
The most intense of these are the lines at 431 and 478 keV from 7 Be* and 
7 Li*. The cross section for 7 Li production in ctct reactions is shown in 
Figure 14 (Kozlovsky and Ramaty 1974b). The cross section for 7 Li* pro- 
duction is about half of the total 7 Li production independent of energy. 
Similarly, the cross section for 7 Be* should also be about half of the 
total 7 Be production; and because 7 Be and 7 Li are produced by mirror 
reactions their production cross sections should be about equal, even 
though no data for 7 Be production in act reactions is available. 

Using these cross sections we find that the intensities of the lines 
at 478 keV and 431 keV are approximately the same as the intensity of the 
4.43 MeV line. Their Doppler widths, however, are about 30 keV (Kozlovsky 
and Ramaty 1974a); therefore, they should not be observed individually 
but rather as a broad spectral feature. Because of this reason and the 
higher continuum emission around 0.5 MeV, it is more difficult to observe 
these lines than the 4.43 MeV line. 

The gamma- ray lines at 0.72, 1.99, 3.62, 3.84, 5.3, 6.33, and 6.7 MeV 
listed in Table 2 are produced by spallation reactions. We have used the 
relative cross sections of Clegg et al. (1961), Foley et al. (1962) and 


Zobel et al. (1965, 1968) 
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to scale the intensities to the spallation-produced contribution of the 
4.43 MeV line at P c = 150 MV assuming similar energy dependences of the 
excitation functions. Similarly, for the lines at 0.845, 1*24, 1.38, 

1.78, 2.75 and 7.12 MeV, which are produced by direct excitation, we 
have used the cross sections of Zobel et al. (1968) and McGowan et al. (1969, 
1970) to scale the intensities to the direct excitation contribution 
of the 4.43 MeV line. 

The 1.63 MeV line results from both direct excitation and spal- 
lation reactions. The cross section for the reactions 30 Ne(p,p ') 3 °Ne* 1 
and 14 N(p,p 1 ) 14 N* a ° 94 are about 0.5 and 0.15 times the l2 c* 4 ° 4 3 direct 
excitation cross section, respectively (Oda et al. I960, Ascuito 1972), 
and the 3.94 MeV level decays 967. of the time to the 2.31 MeV level. The 
cross sections for the spallation reactions on 16 0 and 24 Mg which also 
lead to 1.63 MeV photons are from Zobel et al. (1965, 1968). 

The 2.31 MeV line also results from a variety of processes. The 2.31 
MeV level in 14 N can be populated both by the decay of l4 0 (99.4% of the 
time) and 14 N* 3 * 94 . The cross section for producing the positron 
emitter 14 0 is discussed in Section 7, and the cross section for the 
excitation of the 3,94 MeV level in 14 N was discussed above. The cross 
section for the direct excitation of the 2.31 MeV level is about 5 to 
10% of that of the 3,94 MeV level (Oda et al. 1960), The cross section 
for the reaction 16 0(p, 2pn) 14 N* 3 " 3l is given by Zobel et al. (1968) . 
Combining all the above processes, we find that the contribution of l4 0 
is about 70% of the total. We note that since the half life of 14 0 is 
70*5 sec, the 2.31 MeV line is delayed with respect to the other prompt 


lines of Table 3. 
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In summary. Table 3 shows the most important prompt gamma-ray 
lines that could be observed from solar flares. The most prominent line 
is that at 4.4 MeV, but because Doppler broadening of the 6.14 MeV and 
6.33 MeV lines makes them essentially indistinguishable, a feature at 
^ 6.2 MeV could be as intense as the 4.4 MeV line. Its width, however, is 
about 300 keV, while the Doppler width of the 4.4 MeV line is only ~100 
keV. The lines at 0.43 MeV and 0.48 MeV are also quite intense, but 
because of their larger relative Doppler widths (~3G keV) and higher 
underlying continuum emission, these lines are more difficult to observe. 
The line at 2.31 MeV, even though close to the 2.225 MeV line, makes a 
negligible contribution to line emission at this energy (~3%) . The iron 
lines at 0.85 and 1.24 MeV turn out to be surprisingly strong. These 
lines should be relatively easy to observe with a high resolution detector, 
since their Doppler widths are only a few keV. 
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5. Continuum Emission 

The radiation mechanisms that can produce continuum emission in 
the x-ray and gamma-ray regions are bremsst rah lung, Compton scattering 
and synchrotron radiation but the latter two are not important in solar 
flares. The magnetic field in flares are not expected to be larger than 
about 1000 gauss and the energies of the accelerated electrons do not 
exceed about 100 MeV. Therefore, synchrotron radiation is not expected 
to produce photons of frequencies greater than a few times 10 13 Hz and hence 
this mechanism is unimportant in the x-ray and gamma-ray regions . 

Similarly, Compton scattering (Shklovsky 1965) is also not an import- 
ant radiation mechanism, as can be seen from the following argument. 

In the solar atmosphere, the largest energy density in radiation fields 
is that due to the photospheric emission at T ~ 6000* K. Since the energy 
density in this field close to the photosphere is only about 1 erg cm“ 3 
which is only a fraction of a percent of the energy density in flare 
magnetic fields, any accelerated electron will lose at least 100 
times more energy by synchrotron emission than by Compton scattering. 

The synchrotron energy loss is observable in the microwave region. If the 
energies of the accelerated electrons are smaller than about 100 MeV, the 
Compton scattered photons will have energies below ~ 20 keV; since the 
measured energy flux in these x-rays is much larger than that observed in 
the microwave region, only a very small fraction of this flux can be due to 
Compton scattering of visible photons. Another possible Compton process 
would be the scattering of x-rays below about 20 keV into the gamma-ray 
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region by electrons of a few MeV. But the energy density in these 
x-rays is smaller than in the visible radiation (~ lC^ 2 erg ctir 3 for 
the 1972, August 4 flare, Dere et ai* 1973). Hence also in this case 
the Compton losses are quite negligible, 

We consider now the production of bremss trahlung in solar 
flares. In the nonre lat ivist ic region, most of the radiation is. pro- 
duced by the interaction of accelerated electrons with ambient ions. 

The inverse process, in which accelerated protons interact with ambient 
electrons has also been considered (Boldt and Serlemitsos 1969). If 
the accelerated protons have the same velocity as the accelerated 
electrons, proton-electron bremss trahlung has the same cross section as 
electron-proton bremmstrahlung. But for the one flare (1972, August 4) 
for which we know the proton flux from the gamma-ray observations, it 
turns out that bremss trahlung due to proton-electron interaction is 

negligible in comparison with that due to electron-proton interactions. 

Bremss trahlung can also be produced in electron-electron interactions * 
This process is negligible in the nonre lat ivis tic region but becomes 
comparable to e iectron-proton bremmstrahlung in the relativistic domain 

(Akhiezer and Berestetskii 1965). : 

The differential cross section for electron-proton bremss trahlung 
is given by (Heitler 1954, Koch and Motz 1959) 
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where e and y are the photon energy and initial electron energy, 
respectively, both in units of the electron rest mass energy, 
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is given by an equation similar to equation (1), 

q(e) = n J dE N(E) cp do. (E, € ), (26) 

de 

where E is the kinetic energy of the electron. This equation has 
been evaluated numerically (X. Bai, private communication 1974) 
assuming a power-law kinetic energy spectrum given by equation (14) and the 
results are presented in Figure 15 for various spectral indexes s. This 
figure is for electron-proton bremsst rah lung only. In the subsequent sec- 
tion we shall add the bremsst rah lung due to electron-helium interactions, 
which, for an assumed He/H ratio of 0.07, increases q (E) by a factor of 
1.28. In the nonrelativistic region q (E) can be approximated by a 
power-law proportional to E S °* 6, in the ultrarelativistic region 
q(E) is proportional to fT s+0 ’ 1 . 
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6 . Accelerated Particles at the Sun 

The energy spectrum of the accelerated nuclei at the Sun can be deduced 
by comparing the calculated and observed ratios of the intensities of the 
strongest prompt line at 4.43 MeV to that of the neutron capture line at 
2.225 MeV; and the total number of accelerated nuclei can be obtained by 
comparing the calculated and observed absolute intensities of these lines. 

The energy spectrum and total number of the accelerated electrons can be 
obtained by comparing the observed x- and gamma-ray continuum fluxes with the 
calculated bremsstrahlung intensities. Let us first consider the nuclei. 

The ratio of the intensities of the 4.4 MeV and 2.2 MeV lines is 
given by 



where ( l2 c* 4 ° 43 /n) is the ratio of the yield of the excited nucleus 
i^C* 4 ' 43 to the neutron yield shown in Figure 13, and f is the 2.2 MeV 
photon yield per neutron given in Figure 9. By combining these two figures, 
we can calculate 0^ 4 /^2 2 from equation (27). The results are plotted 
in Figures 16 and 17, for 3 He/H = 0 and 3 He/H = 5xl0" 5 , respectively. The 
results of Figure 16 differ only slightly from those of Figure 10 of Ramaty 
and Lingenfelter (1975a, b) who have used a constant value for the 2.2 MeV 
photon yield per neutron (f = 0.2) instead of the model-dependent results 
of Figure 9. 

Let us compare now the calculations with the data. The observed 
(Chupp et al. 1975) ^ 4 , 4 / ^2.2 ratio of 0.11 + 0,04 for the 1972., August 4 
flare is also shown in Figures 16 and 17. As can be seen, assuming power- 
law spectra for the particles in the flare region, their spectral index s 
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should lie between the values of 1.9 + 0.2 deduced for the thin-target model 
and 3 + 0.3 for the thick-target model if there is no 3 He in the photosphere; 
or between 1.7 +0.2 and 2.7 + 0.3 for these models if the photospheric 
3 He/H ratio is 5x10" 5 . Similarly, assuming exponential spectra, the implied 
P 0 's should lie between 110 + 30 MV for the thick-target model and 170 + 50 MV 
for the thin-target model if 3 He/H = 0; or between 165 £ 55 MV and 250 + 80 MV 
for these models if 3 He/H = 5xl0" 5 . These spectral parameters are summarized 
in Table 4. Note, however, that these implied spectra become somewhat 
steeper if we use a larger He abundance, since, as discussed in Section 2, 
this increases the neutron production without changing the 12 C* 4 * 43 

Knowing the spectral index for the various models, we can now deduce the 
total number of protons at the Sun. In Figures 18 and 19 we show the instant- 
aneous and total 2.2 MeV fluxes for these models, for a He/H = 0 and 
3 He/H = 5xl0~ 5 , respectively. These figures are obtained from equations (21) 
and (22), where q, Q and f are from Figures 4 and 9. 

In the thick-target model the time -integrated photon flux from the flare 
determines the total number of accelerated particles that interact and stop 
at the Sun. According to Chupp et al. (1975), the 2.2 MeV intensity for the 
1972, August 4 flare was about 0.3 photons cm -2 s -1 . This flux is the average 
over the time interval of observation (0623 to 0633 UT) which gives a total 
flux of about 180 photons cm -2 . Since the detector on 0S0-7 was eclipsed 
by the Earth before the termination of the gamma-ray event, this is a lower 
limit to the total flux from the flare. If we assume a duration of ~10 3 
seconds for the event, as indicated by the hard x-ray data (van Beek et al. 
1973), then the total flux of 2.2 MeV photons was ~300 cirT 2 . Then from 
Figure 18 we find that N p (>30 MeV) = 8xl(P 2 for s = 3, and N p (>30 MeV) 
= 1.0 x 10 33 for P Q = 110 MV. The corresponding values for 3 He/H = 5xl0“ 5 
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are about the same, because the effect of a lower f is approximately cancelled 
by a higher neutron yield* Therefore, the total number of protons above 
30 MeV released downward into the Sun in the thick-target model is about 
10 33 independent of the spectral form or the 3 He/H ratio. 

For the thin- target model, the instantaneous gamma-ray observations 
determine the product of the ambient proton or hydrogen density, n H , and 
the instantaneous number of the accelerated particles in the interaction 
region. From Figure 18 we obtain that n H N p (>30MeV) 21 SxlO^cnT 3 
for s = 1,9, and n H N p (>30MeV) ^ 8xl0 43 cm’ 3 for P Q = 170MV„ As before, 
these values are not significantly effected by the photospheric 3 He/H ratio* 

We consider now the electrons* The observed x and gamma-ray continuum 
flux for the 1972, August 4 flare is shown in Figure 20, The lower and 
upper envelopes of the shaded region are the best fits to the observed 
x-ray flux at about 0626UT and 0630UT, respectively (van Seek et al. 1973), 

The data points are from Suri et aL (1975) and they represent the average 
gamma-ray flux over the time interval 0623UT to 0633UT, Thus, even though 
the x-ray and gamma-ray spectra were observed by different detectors on 
different satellites, the observation times were about the same, and there is 
reasonable agreement between the different observations* 

Suri et al* (1975) have discussed some of the implications of the 
photon spectrum shown in Figure 20 „ The change of slope around 100 keV has 
been observed previously in other solar flares (Frost 1969), and it is probably 
due to a break in the electron spectrum caused by the acceleration mechanism 
that operates in this energy region, A more novel feature of the spectrum in 
Figure 20 is the flattening at ~500 keV which leads in the MeV region to an 
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excess photon flux over the simple extrapolation of the observed flux at 
lower energies c It is worth noting that a similar flattening was probably 
observed for the 1967, May 23 flare by Gruber et al. (1973) * The absolute 
photon flux in the MeV region for this flare was almost the same as for the 
1972, August 4 event c 

Suri et al. (1975) have suggested that this excess could be due to a 
separate population of electrons which would be accelerated by a second stage 
of acceleration 0 Independent of the origin of these electrons, we proceed 
now to calculate their number and spectrum from the observed gamma-ray 
continuum. 

The line A in Figure 20 represents the instantaneous photon flux 
obtained from the electron spectrum A in Figure 21 by using equation (26) 
and the bremsstrahlung cross section given by equation (24) □ Similarly, 
the lines B and C represent the photon fluxes obtained from the electron 
spectra B and C 0 The spectrum A has a break at about 140 keV which produces 
the observed break in the photon spectrum at about 100 keV and it has a high 
energy cutoff at 1 MeVo Both the electron spectra B and C can produce the 
observed gamma -ray continuum above ^0.5 MeV c 

The dashed lines in Figure 21 are some representative proton spectra 
in the thin- target model obtained by using the parameters of Table 4, 

Because the gamma-ray lines which determine these parameters are quite 
insensitive to the proton spectrum below seyeral MeV, while the gamma-ray 
continuum which determines the electron spectrum has been measured only up 
to several MeV, there is no common energy region where both the proton and 
electron spectra are well determined* Nonetheless, the results of Figure 21 
seem to suggest that the proton- to-electron ratio at the same energy ranges 
from about 10 to 10 2 . The analytic forms of the spectra shown in Figure 21 
are given in Table 5. 
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Comparison of these implied particle spectra in the solar flare region 
with the proton spectrum observed in the interplanetary medium from the 
1972, August 4 flare is complicated by the possibility that the latter 
spectrum may have been significantly modified by acceleration in the inter- 
planetary medium. 

To obtain an estimate of the number of protons released into the inter- 
planetary medium, the information obtained from the gamma rays in the thin- 
target model can be combined with data on the path length traversed by the 
nuclei before their escape from the interaction region. Such information 
can be obtained from studies of deuterons and helium-3 abundances in accel- 
erated nuclei from flares. The 3 H and 3 He observations from the 1972, 

August events (Webber et al. 1975), when compared with calculations of the 
production of these isotopes in nuclear reactions (Ramaty and Kozlovsky 1974) 
imply that the amount of matter traversed by relativistic particles is about 
1-5 g cm“ 3 . This means that the product n^t t is about 3xl0 13 cnT 3 s, where tj 
is the interaction time of the particles at the Sun. If is also interpreted 
as the escape time of the particles from the interaction region, then 
from the values of determined above, we find that the protons were 

released into the interplanetary medium at an average rate N^C^OMeV)/^ 
varying from about 1.6 to 3X30 30 protons s’" 1 . The total number of protons 
released is the product of this rate and the acceleration time T. As indicated 
by the x-ray data (van Beek et al. 1973), T is about 10 3 seconds; therefore 
the total number of protons released above 30 MeV is between about 1.6 to 
3x10 s3 . These numbers are larger by only about a factor of 2 to 3 than the 
number of protons released downward into the Sun in the thick- target model. 

The number of protons released from the flare should be compared with 
estimates of the total number of protons in the interplanetary medium as 
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obtained from charged particle observations. Ramaty and Lingenfelter (1973a) 
have estimated this quantity by taking the observed density of protons near 
Earth and multiplying it by a storage volume in the interplanetary 
medium which they took to be ~ 10 39 ciri 5 . Using this volume and the peak 
proton flux greater than 30 MeV measured by Kohl et al. (1973), we get about 
IQ 34 protons « This number is larger by 2 factor of 3 to 6 then, our estimate 
for the number of protons released in the thin-target model * But since 
the measured peak proton flux could consist to a large degree of particles 
accelerated by shocks in the interplanetary medium, it appears that there 
is no real discrepancy between the number of protons released for the Sun 
as deduced from the gamma rays and the number observed in interplanetary 
space . 

We now turn briefly to the processes of energy loss by the accelerated 
particles in the flare region and the problem of whether such processes are 
the principal source of energy for the heating of a flare . Giovanelli (1946) 
first suggested that accelerated electrons caused the ionization and excit- 
ation responsible for the optical line emission, while Gordon (1954) proposed 
that accelerated protons were primarily responsible. Detailed studies of the 
ionization losses and resulting optical emission processes have been made 
for both accelerated electron (e.g. Neupert 1968; Brown 1973; 

Somov and Syrovatskii 1974) and proton (e 0 g 0 Dubov 1963; Najita and Orrall 
1970; Svestka 1970) interactions with the ambient gas in the flare region. 

Although estimates (Neupert 1968; Syrovatskii and Shmeleva 1972; 

Brown 1973) of the electron energy loss, based on measurements of the 
x-ray emission suggest that sufficient energy for the optical emission may 
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be available from the electrons alone, an estimate of the relative importance 
of electron and proton ionization losses has not previously been possible 
because of the lack of information on the number and spectrum of accelerated 
protons in the flare region 0 The gamma ray measurements of Chupp et al* 

(1973) however, now allow us to at least place some limits on the ratio of 
total energy loss rates of accelerated protons and electrons in the 1972, 
August 4 flare o 

For the thin-target case the instantaneous energy loss rate for 
electrons and protons is 

W = J J* dE N (E) cpdE/dx (28) 

where the electron collisional energy loss rate dE/dx for electrons is taken 
from Berger and Seltzer (1964) and for protons from Barkas and Berger (1964) . 

For the thick-target case where essentially all the entire kinetic energy 
of the particle is dissipated in ionization losses the total energy loss 
is simply 

W = J dE N (E) E ( 29 ) 

for both electrons and protons. The energy dissipated by electrons in 
bremsstrahlung and by protons in inelastic collisions is neglible fraction 
of the total kinetic energy for the flare accelerated particle. 

The minimum energy loss by accelerated electrons in the 1972, August 4 
flare can be directly determined from the x-ray emission measured by van Beek 
et al . (1973) down to about 20 keV , shown in Figure 20. Below this energy 
thermal emission from the high temperature region of the flare masks the 
nonthermal emission. As discussed above, assuming that the emission >20 keV 
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results from accelerated electron brems st rah lung ^ the required total electron 
energy spectrum for the thin-target case is found to be 

N(E) = 2.1xl0 44 E(MeV)' 2 « 5 /n, (30) 

in the energy range from ~20 to ^140 keV. The greatest contribution to the 
energy loss comes from these lowest energy electrons. Using this spectrum 
in equation (28) we find an instantaneous electron energy loss rate 
W (>20 keV) - 3X10 28 ergs/sec. Also, as discussed above, if we assume a time 
of the order of ~10 3 sec for the duration of the flare we find a total electron 
energy loss W (>20 keV) fsOxIO 31 ergs. A comparable total energy loss is 
also required to account for the observed emission by electron bremsstrahlung 
in the thick-target case. As in other flares the energy loss from these 
electrons exceeds the total radiation energy of the flare, comprising an 
estimated (Zirin and Tanaka, 1973) total optical emission energy in of 
~2xlO ao ergs, with comparable energy emitted in Lyman a and in high temperature 
thermal emission in soft x-rays. A large fraction of heating energy, however, 
may remain in the plasma cloud blown off in the explosive phase of the flare. 

The most direct estimate of the energy loss of accelerated protons in 
the flare can be made from the prompt gamma-ray line emission from 13 C* 4 ‘ 4 . 

In the thin-target case, the ratio of the energy loss rate to the 4.4 MeV 
line intensity at Earth is obtained by dividing equation (28) by equation 
(12) times 4 ttr 3 , where R = 1A.U. In equation (28) for power-law spectra, 
we have assumed that below an energy E c , N(E) is constant, i.e. <N(E) = 
kE for E>E C , and N(E) = kE c for E<E C . No such cutoff is assumed for 

o 

exponential spectra. The ratio Wp /0 4 , 4 shown in Figure 22 is independent 
of the ambient gas density; it depends only on the spectral index, and for 
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power-law spectra, on the cutoff energy E c « As can be seen from Figure 10, 
the production of ^C* 4 * 43 is independent of the assumed cutoff energy if 
E c is less than a few MeV. 

For power- law spectra in the range 1,7 < s < 2*1 determined 

above from the ratio 2,2 MeV to 4 0 4 MeV line intensities, we see from 
Figure 22 that ft/0 4#4 for E c = 2„5 ranges from l.AxlO 28 to Ic^xlO 38 ergs 
cvP /photon* Thus the observed 0 4 4 of ~3xl0" 2 photons /(tra 2 sec) requires a 
minimum proton-energy loss rate of ^(5 + l)xl0 26 ergs/sec . There is no 
reason to believe, however, that the proton spectrum would cut off just below 
the l2 C* 4 ° 4 threshold. If the same spectral index held down to 0.5 MeV or 
0.1 MeV, for example, the required proton energy loss rate, derived from 
Figure 22, would then be (3 + 1.5)xlO s7 or (2 + l o 6)x!0 28 ergs/sec respect- 
ively. The latter value is comparable to the minimum electron energy loss 
rate of 3xl0 28 erg/sec* For the exponential spectra, on the other hand, the 
range 120 MV < P Q < 300 MV gives a proton energy loss rate of only (4 + lJxlO 26 
erg/sec independent of any assumed cutoff energy. Thus we see that the 
uncertainty in the low-energy cutoff of the spectrum does not allow us to 
determine whether or not proton ionization losses make a significant contri- 
bution to the heating of the flare region. 

\ 

A similar situation prevails for the thick-target case. The ratio 
of the total proton energy loss to time -integrated 4*4 MeV gamma-ray line 
emission, W/0^ 4 , obtained by dividing equation (29) by 4 TTr 3 times equation 
(13), is also shown in Figure 22. Here, for the thick-target power-law 
spectra with 2.4 < s < 3.2, we find a total proton energy loss 
W & (2 + l)xl0 3 ° ergs for the minimum cutoff E c - 


2*5 MeV or energy 
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losses of ~(1.1 + 0„9)xl0 31 ergs for E c = 0.5 MeV and (7 + 6„5)xl0 31 ergs 
for E c = 0.1 MeV. These are again comparable to the minimum electron loss 
of ^xlO 31 ergs, but the uncertainty in E c permits no more quantitative 
comparison. For exponential spectra, however, with 80 MV < P Q <• 220 MV, 
the total proton energy loss must be only ~(2.1 + 0.6) x 10 29 ergs. 
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7 . The Nature of the Positron Annihilation Radiation 

Positrons can be produced in solar flares from the decay of j f*" mesons 
and radioactive nuclei which result from nuclear reactions of accelerated 
protons and nuclei with the ambient medium. Some solar positrons could 
also be created by accelerated electrons in e + -e" pairs. Annihilation of 
positrons can produce a gamma -ray line at 0.51 MeV . This line was observed 
by Chupp et al. (1973) for both the 1972, August 4 and August 7 flares. 

The most important positron emitters, their half lives and maximum 
positron energies are listed in Table 6 , together with their main pro- 
duction reactions and their threshold energies. The cross sections for 
the production of jr - mesons were given by Lingenfelter and Ramaty (1967). 

The cross sections for producing the radioactive nuclei are shown 
in Figure 23. The cross sections for forming 10 C and i:L C by spallation of 
16 0, 14 N, and 12 C were summarized by Audouze et al. (1967), Meneguzzi et al. 

(1971) and Mitler (1972). To these cross sections we added the low-threshold 
reaction 14 N(p, ct ) 11 C measured recently by Jacobs et a l. (1974). The croS s section 
for producing l2 N by 12 C(p,n) 12 N was taken from Rimmer and Fisher (1968). 

The 13 N producing-cross section is the same as in Lingenfelter and Ramaty 
(1967) except for the low-energy part of the reaction 16 0£p, ) 13 N which 
is from Albouy et al. (1962). The cross sections for the reactions 
14 N(p,n) 14 0 and ls O(p, ) l4 0 are from Audouze et al. (1967) and Meneguzzi 
et al. (1971). The cross section for 15 0 production in proton induced 
reactions is the same as in Lingenfelter and Ramaty (1967). In the numerical 
calculations we have added the low-threshold reaction 12 C(o,,n) l5 0 which has 
been discussed in the neutron-production section. The cross section for the 
reaction 16 0(cv,n) l9 Ne has not been measured. As discussed in Section 3, we 
assume that this cross section is the same as that of the reaction 12 C(o',n) 15 0 
which was measured by Nelson et al. (1963). 
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Using these cross sections we calculate from equations (12) and (13) 
the positron emitter production for the thin and thick-target models. 

By way of illustration the instantaneous production rates of positron- 
emitters in the thin-target model for power-law spectra are shown in 
Figure 24. As can be seen, for flatter spectra the principal positron 
source is tt* mesons, but for steeper spectra 11 C is the largest source and 
the pion contribution is negligible. This effect can also be seen in 
Figure 25, where the ratio of the pion yield to the total positron-emitter 
yield (including both pions and radioactive nuclei) is plotted as a function 
of s or P Q for both interaction models. We see that, as for neutron 
production, the efficiency of T x + production is larger in the thick-target 
model than in the thin-target model. The ratios of the total yield of 
positron-emitters to the neutron yield are shown in Figure 26 for the 
various interaction models. 

The intensity of the 0.51 MeV line, however, depends not only on the 
number of positron emitters produced, but also on the decay rate of the 
positron emitters and on the annihilation rate of the positrons. The half 
lives against radioactive decay of the various positron emitters are shown 
in Table 6. Thus, in the early phases of particle acceleration and inter- 
action on time scales less than 20 minutes, is a less important source 

than the shorter lived (1 to 2 minutes) positron emitters l4 0 and 15 0. At 
later times, however, 11 C can provide delayed positrons after significant 
nuclear interactions have ceased. 

In addition, the annihilation rate of the positrons depends on the 
density, temperature and state of ionization of the ambient medium. As 
discussed previously (Stecker 1969, Ramaty and Lingenfelter , 1973b) in a 
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low-density neutral medium, positrons nearly always annihilate via positronium 
formation. Only if the ambient density exceeds about 10 15 cm -3 do col- 
lisions destroy the 3 S state of positronium at such a rate that free 
annihilation becomes more important than positronium annihilation (leventhal 
1973). Positron annihilation from a bound state of positronium results 
in an asymmetric 0.51 MeV line, since 75% of the time positronium annihilates 
from the 3 S state into 3 photons of energies less than 0.51 MeV, instead 
of 2 photons at precisely this energy. The 0.51 MeV line from the Sun, 
however, appears to be symmetric. 

As we noted, three-photon annihilation could be reduced by collisions, 
but the required density of ^ 10 15 cm -3 is quite large. However, if the 
ambient medium is ionized, the rate of positronium formation is greatly 
reduced. In Figure 27 we show the rate of positronium formation and free 
annihilation per positron in a hydrogen plasma of unit density as a function 
of its temperature (C. Werntz and C.Crannell, private communication 1973). 

As can be seen, if the temperature is greater than about 1C)P°K, most 
of the positrons annihilate without forming positronium. This is not an 
unreasonable temperature for the annihilation region. From the observed 
upper limit on the width of the 0.51 MeV line, however, the temperature in 
the annihilation region should be less than ~ 10 7o K (Chupp et al. 1975). 

According to the observations of the 1972, August 4 flare, the rise 
time of the 0.51 MeV line was, within instrumental errors, similar to or 
perhaps even shorter than the pise time of the 2.2 MeV line (Chupp et al. 
1975). The latter was about 100 seconds, consistent with the expected lag 
between the production of the neutrons and the formation of the 2.2 MeV 
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line (see Table 2). Assuming that the time dependence of the production 
rate of the positron emitters is the same as that of the neutrons, it 
follows that the lag between the production of the positron emitters and 
the formation of the 0.51 MeV line should also not be longer than about 
100 seconds. From Figure 27, this result implies that the density of the 
ambient medium in the annihilation region is at least 10 12 cnT 3 . Further- 
more if the positrons are of nuclear origin, then from Table 6 it follows 
that they should mainly result from tt + mesons and short lived radioactive 
nuclei ( 15 0, 14 0). 

The observed average flux in the 0.51 MeV line for the flare of 1972, 
August 4 was about 0.06 photons cnT 2 s- 1 (Chupp et al. 1975). Hence the 
observed 51^2. 2 ratio was about 0.2. 

Let f' be the 0.51 MeV photon yield per positron defined in the same 
way as the 2.2 MeV photon yield per neutron, f, in equation (21). As 
discussed in Section 3, f ranges from about 0.1 to 0.2. The maximum value 
of f 1 is 2. But because part of the positrons can escape from the Sun 
before they annihilate, and a fraction of the positrons can be trapped at 
the Sun in low density regions where the annihilation time is long, f can 
be considerably less than 2. 

According to our discussion in Section 6, for power law spectra s~2 in 
the thin-target model, and s~3 in the thick-target model. For these spectral 
parameters, from Figure 25 we see that the bulk of the positron emitters 
are TT i mesons, and from Figure 26 we get that the positron emitter-to- 
neutron ratio is about 0.2. The observed ratio, 0 q 5x /^2 2 — 0.2 then 

implies that f' ~ f, i.e. the 0.51 MeV yield per positron is about the same 
as the 2.2 MeV yield per neutron. In view of the uncertainties involved 
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in the deductions of f', this is not an unreasonable result. We can also 
deduce f* for the various values of P 0 obtained in Section 6, and in all 
cases we find acceptable values (0.1 < f' < 2). For these exponential 
spectra the contribution of tt + mesons to the total positron production is 
greater than about 50%, in all cases, except in the thick-target model with 
no 3 He in the photosphere where they contribute only about 25% of the 
positrons. A definite test of the possibility that the bulk of the positrons 
are due to rr + mesons would be the observation of gamma rays from tt° decay. 

We give estimates of the expected flux of such gamma rays in Section 8, 

As mentioned above, prompt 0.51 MeV photons could also result from 
positrons produced in e + -e~ pairs. The cross section for pair production 
by relativistic electrons on ambient protons is given by (Heitler 1954) 

Q pair ~ 2^(* ro ) 2 [ ln ( H ] 3 < 31 

Using this cross section, the instantaneous pair production rate was 
calculated (T. Bai, private communication 1974) from equation (12) for the 
thin-target model with power-law spectra. The results are shown in Figure 
28 for various spectral indexes s and high-energy cutoffs, E T , of the assumed 
electron spectrum. 

As we have seen in Section 6, the continuum gamma-ray emission from 
the flare of 1972, August 4 is consistent with an electron spectrum which 
above about 1 MeV is given by curves B or C in Figure 20. If we consider 
the spectrum B, and assume that the 0,51 MeV photon yield per positron, f', 
has its maximum value of 2, then from Figure 28 we find that for the 1972, 
August 4 flare the 0.51 MeV photon flux from pair production ranges from 
about 8x10" 5 to 6X10 -4 photons cm -2 s“ 1 , depending on the assumed high 



energy cutoff* These values are smaller by factors of 10 3 to 10 3 than the 
observed flux of about 0 o 06 photons cnT^s" 1 (Chupp et al* 1975)* Pair 
production therefore accounts for less than 1% of the observed 0*51 PleV 
line intensity* 



High Energy Gamma Rays and Neutrons 


In addition to line emissions, nuclear reactions in solar flares also 
produce TT ° mesons which decay into high energy gamma rays P Each tt° meson 


yields two photons with energy spectrum centered about 67,5 MeV. Using the 
tt° production cross sections in proton-proton (Lingenfelter and Ramaty 
1967) and proton-^ particle reactions (Stecker 1970), we have calculated 
the yields of TT° mesons in the thin and thick-target models for power-law 
and exponential spectra. The results are shown in Figure 29« The total 
TT°-decay gamma-ray flux can be obtained from this figure by using the formula 


0 


TT 


= 0 


- (^) 


7 -\ 


(32) 


where 0 p o;3 is the observed 2.2 MeV line intensity and f is the average photon 
yield per neutron given in Figure 9. 

No high-energy gamma rays have been observed from solar flares. Nonetheless, 
the following values of 0^° can be calculated for the 1972, August 4 flare: 

For 3 He/H = 5xlCT 5 , the parameters of Table 4 yield: 


3xlCr3 




0.13 

for 

the 

thick-target 

exponential model. 

0.035 

£ 

0 o 

TT 

£ 

0.37 

for 

the 

thin-target 

exponential model, 

0,13 


0 TT° 


0.42 

for 

the 

thick-target 

power-law model, and 

0.25 

£ 

0 TT° 

£ 

0.7 

for 

the 

thin-target 

power-law model * 


These fluxes should be compared with the bremsstrahlung that is expected at 
high energies from the electron spectra shown in Figure 21, For the spectra 
B and C we get 0(>3O MeV) ~ 0-07 ph ctrT 3 s -1 and 0(>3O MeV) ~ 10” 3 ph cm -2 s" 1 , 
respectively. We see that high-energy gamma rays from 110 decay could be 
observable from solar flares. Moreover, the characteristic broad maximum 
of these gamma-rays' energy spectrum makes them easily distinguishable from 
a possible bremsstrahlung continuum. 
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Finally, we also present the expected-high energy neutron fluxes at 
Earth for our various interaction models. In Figure 30 we give the ratio of 
the time -averaged neutron flux 0^ to the time -averaged 2.2 MeV line intensity 
for these models. The time dependence of these two fluxes, however, are 
quite different, as was discussed by Lingenfelter and Ramaty (1967). 
Nonetheless, we see that for certain parameters, 0 n can become as large or 
even larger than 0 g , and hence potentially detectable for large solar 
flares. We note that 0 n in Figure 30 has been calculated under the assumption 
that the neutrons are produced isotropically at the Sun. However, 0 n is 
greatly reduced if the neutrons are produced with initial velocities pre- 
dominantly toward the Sun. 

The detection of gamma-rays from tt° decay or high energy neutrons would 
provide spectral information on accelerated particles in the several hundreds 
MeV /nucleon range. In contrast, the nuclear deexcitation lines and the 2.2 
MeV line which is produced by lower energy neutrons, provide information 
on the primary particles mainly in the 10 to 100 MeV/nucleon region. 
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9 a Summary 

The observed gamma-ray lines from the 1972, August 4 flare, at 0.5, 

2.2, 4.4 and 6.2 MeV are due to positron annihilation, neutron capture on 
hydrogen, and deexcitation of excited states in 13 C, 16 0 and 1S N, respec- 
tively. The strongest line is at 2.2 MeV. It is due to fast neutrons 
produced by nuclear reactions of flare accelerated particles with the ambient 
solar atmosphere. These neutrons are thermalized and captured by ambient 
protons in the photosphere to produce deuterons and 2.2 MeV gamma rays. 
Photospheric 3 He competes with the protons in capturing neutrons 0 Because 
captures on 3 He do not lead to photon emission, the observation of 2.2 MeV 
line emission from the Sun implies that the 3 He abundance in the photosphere 
cannot be much larger than that observed in the solar wind ( 3 He /H ^ 5xlCT s ). 

We have evaluated in detail the yield of neutrons and excited nuclei 
from nuclear reactions of accelerated particles with the ambient solar 
atmosphere. For the 1972, August 4 flare the neutrons are produced mainly 
in pcv and <yp reactions by accelerated particles with energies greater than about 
30 MeV/nucleon 0 The observed gamma rays at 4.4 and ^6.2 MeV are principally 
due to proton induced interactions. Reactions induced by fast nuclei lead 
to Doppler broadened lines which cannot be distinguished from the continuum. 

We have determined the relative importance of other prompt lines that might 
be observed from solar flares in addition to the 4.4 MeV line and the 
6.2 MeV feature. The most promising candidates are at 0.85, 1.24, 1.63, 1.78, 
2.31 and 5.3 MeV. Of particular interest are the iron lines at 0.85 and 1.24 
MeV which have very narrow Doppler widths of only a few keV. 

The positrons which produce the 0 o 5l MeV line are due to TT+ mesons and 
radioactive nuclei. We have evaluated in detail the production of mesons 



- 46 


and radioactive nuclei in nuclear reactions of accelerated charged particles 
with ambient nuclei. The relatively prompt nature of the 0.51 MeV line 
seems to favor positron production from tt + decay in the initial phase of the 
1972, August 4 event. The 0.51 MeV line emission observed late in the 
1972, August 7 event is very likely due to delayed positrons from radioactive 
nuclei. The spectral region around 0.5 MeV is c OTTin 1 ic3 1 0 d by the possibility 
that some of the positrons decay into 3 photons from the triplet state 
of positronium, and that a strong feature at ~.46 MeV is produce by the 
7 Li and 7 Be lines. A detector with very good energy resolution is required 
to resolve the various structures in this energy region. 

From the comparison of the calculated and observed ratios of the line 
at 4.4 MeV to the line at 2.2 MeV it is possible to deduce the spectrum 
of the accelerated particles in the flare region. For the 1972, August 4 
flare the spectral index s defined in equation (14) is 1.9 + 0.2 for the 
thin- target model and 3+0.3 for the thick- target model if there is no 
3 He in the photosphere; or 1.7 + 0.2 and 2.7 + 0.3 for these models if the 
photospheric 3 He/H ratio is 5x10" 3 . The characteristic rigidity P Q defined 
in equation (15) is 170 + 50 MV for the thin-target model and 110 + 30 MV 
for the thick-target model if 3 He/H = 0; or 250 + 80 MV and 165 + 55 MV for 
these models if 3 He/H = 5xl0~ 5 „ The total number of protons above 30 MeV 
released downward into the Sun in the thick-target model is about 10 33 . For 
the thin-target model about 2x10 s3 protons of energies greater than 30 MeV 

escape from the flare region if relativistic particles traverse about 
^ 8 CID 3 as deduced from deuteron and helium— 3 abundances in the accel- 

erated particle fluxes. The energy loss rate of protons with energies 
greater than 2.5 MeV is about SxlO 26 erg s' 1 and the total energy 



deposited by these protons at the Sun is ^SxlO 29 ergo However, if the proton 
spectrum as deduced from the gamma-ray observations is extrapolated down to 
0.1 "MeV, the protons deposit an energy comparable to that deposited by 
electrons of energies greater than 20 keVo From the x-ray observations, 
the rate of energy loss of these electrons is ^3xl0^ 8 erg s~ 1 , and the total 
energy deposited by them is ^xlO 31 erg. 

For the flare of 1972, August 4, a gamma-ray continuum has been observed 
at energies greater than about 0*5 MeV, and it appears that it is produced 
by a separate population of electrons in this energy range* The number of 
these electrons is lower than the number of protons of the same energy by 
about a factor of 10 to 10 2 , but their spectrum could be the same as that 
of the protons. This result may have significant implications on particle 
acceleration in solar flares. 

Finally, high energy neutrons and gamma -rays from tt° decay could be 
detected from flares such as the 1972, August 4 event. We expect a 
minimum TT° -decay gamma-ray flux of 3xl0“ 3 photons ciir 2 s" 1 , and, if the 
neutrons are produced isotropically at the Sun-, a time integrated high 
energy neutron flux comparable to the time-integrated 2*2 MeV line intensity. 
The detection of these high-energy gamma -rays -and neutrons would provide 
spectral information on the accelerated particles in the region of several 
hundred MeV. The presently observed nuclear excitation lines and the 2*2 
MeV line provide information on particles mainly below 100 MeV* 

In closing, we wish to note that essentially all the observational data 
used in this paper has been obtained from one large solar flare* As we have 
tried to show, this very limited data already contains a great deal of 
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information not only on the physics of nuclear interactions in flares but 
also on solar physics in general „ We hope, therefore, that more data will 
be forthcoming during the next solar maximum, and that solar gamma-ray 
astronomy will be a vigorously explored discipline of solar physics. 
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TABLE 1 

NEUTRON PRODUCING REACTIONS 


Reaction Threshold (MeV/nucleon) 


1. p-PH -> 

n-fp+'n* 

292.3 

2. p+^Ue -> 

^He+p+ri+On) 

25.7 

-*• 

2 H*f2 pfn+ (jt ) 

32.6 


3p+2n+( TT ) 

35.4 

3. p-P 2 C 

n+ • • o 

19.6 

p-P 3 C - 

n - ! - * o o 

3.2 

4. p-P 4 N -> 

rH- • o • 

6.3 

5 . p-P 6 0 -» 

11+ a t o 

16.6 

p-P 8 0 -*■ 

n+ a o • 

2.5 

6. p-P°Ne -*■ 

n - ^ o • • 

15.9 

7 . p-P 6 Fe 

n-K o • 

5.5 

8. a + 4 He -> 

7 Be+n 

9.5 

9. a-P 2 c - 

n+ «to 

ro 

© 

oo 

0..+13C -► 

XI t g O • 

— 

10.- a-P 4 N - 

n+ 009 

1.5 

11. a-P 6 0 -* 

n+ o o • 

3.8 

a-P 8 o -> 

rH - o o o 

0.21 

12. a-P°Ne -» 

IL’i" » * » 

2.16 

a-p2 Ne _* 

n-K o • 

0.15 

13. a -P 6 Fe -► 

n-f . * , 

1.37 

14. cy -P s Mg -* 

n+« ♦ 4 

___ ’ 

15. a -P 6 Mg - 

n+. , , 



16. a + ?9 Si -» 

. * 

n a*} 


exoergic 



TABLE 2 

Most probable neutron capture densities, capture times, 
the neutron mean life. . 


E n (MeV) 

, -3 n 

<n> (cm ) 

He 3 /H = 5xl0' 5 

1 

7x10 16 

119 

10 

1 . 2x10 17 

69 

100 

3xl0 17 

28 


t , and X = t + Tj , where t. is 
c c d d 

c) X _1 (sec) 

He 3 /H - 0 He 3 /H = 5x10 " 5 He 3 /H = 0 

214 105 173 

125 64 110 


50 


27 


47 



REPRODUCIBILITY OF THE 
ORIGINAL PAGE IS POOR 
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TABLE 3 


Photon Energy 
(MeV) 

PROMPT 

Origin 

GAMMA RAY LINES 

Production Modes 

Approximate 

Re la t ive Intensity 

0.431 

7 Be * 0, 431 g.s. 

4 He ( a ,n) 7 Be* 0431 

1 

0.478 

7 Li*°« 478 -+ g.s. 

4 He( a ,p) 7 Li*°‘ 47e 
4 He ( a , n) 7 Be 127 7 Li*° ♦ 47s 

— y 

1 

0.72 

10 B*° ° 7s -*■ g . s . 

l2 C(p,2pn) l0 B*° * 72 
l6 0(p, ) lo B*°* 72 

0.1 

0.845 

56 Fe *u .84 5 ^ g<s> 

56 Fe (p ,p ' ) eS Fe* 0 * 845 

0.2 


1.24 SS Fe * 2 . 0 S _ S 6 Fe *0 .845 56 Fe ( p f p » )S 6 Fe *2 .0 8 0.2 

1.38 24 Mg*i ° 3? -*• g.s. 24 Mg (p ,p ' ) 24 Mg * 1 ,37 0.1 


1.63 ao Ne* 1,63 -*■ g.s. 20 Ne(p,p') 20 Ne* 1 ° 6a 

i 4 N* 3 * 94 -» N* 2 * 31 14 N(p, p ' ) 14 N*® * 94 

16 0 (p j 2pn) 14 N* 3 ° 94 0.2 

23 Na* 2 *° 7 -> 3 3 Na *0*44 3'4 Mg(pj2p) 23 Na *2.6 7 


1.78 

28 Si *l .78 

g.s. 

2 ?Si(p,p 1 ) 28 Si* 1 * 78 

0.2 

1.99 

11c* 1 099 

-► goS. 

I2 C(p,pn) 11 C* 1 *" 

0.07 

2.31 

1 4 * 3 1 

g*S. 

14 N(p,p , ) 14 N* 3 - 94 -> 14 N* 2 -3l 
14 N(p,n) l4 0* 14 N* 2 * 31 
16 0(p } 2pn) l4 N* a “ 31 

0.3 

2.75 

1 6 q^V 8 • 8 8 

is q*6 .14 

l6 0(p,p') 16 0* 8 * 88 

0.07 

~3 .62 

13C*3 c 68 
6 Li* 3 ° s6 

-* g.s. 

g*S. 

16 0 (p ,3pn) l3 C* 3 068 
12 C(p, ^Li* 3 ’ 56 

0.07 

3.84 

1 3 QVc3 *66 

“-*■ g O S q 

l6 0(p,3pn) 1,3 C* s ° 84 

0.07 

4.43 

l^C* 4 * 43 

g-So 

12 C(p,p') 12 C* 4 * 43 
l3 c(a,a') l2 c* 4 » 43 
l6 0(p,p a ) l2 C* 4 - 43 

1 

~5.3 

15 9*6 .36 
1 5^5 *3e 
1 S |^* 5 *3 1 

gaS. 
g.So 
^ g.So 

16 0 (p,pn) ls 0* 5 “ 3S 
16 0(p,2p) l5 N* 6 * 28 
ls O(p,2p) l5 N* s * 31 

0.3 

6.14 

,■ l 6 o* 6 ° 14 

*+ g«S. 

16 0(p,p , ) l6 0* 6 * l4 

16 0(o.»a') 16 o* 6 * 14 

16 0(p,2p) ls N* 6 ' 3 

0.5 

6.33 

15N * 6 °3 3 

g.s® 

0.5 

~6.7 

11 B *6 «76 

-+ g.s. 

lS C(p,2p) 11 B* 6 * 76 
lS C(p,2p) 11 B* 6 * 81 
l2 C(p,pn) 11 C* 6 * 60 

0.07 

7.12 

1 6 Q * 7 0 1 2 

g.s. 

16 0(p,p') 16 0* 7 ->12 

0.2 
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TABLE 4 


Deduced Spectral Parameters for the 1972, August 4 Flare 


3 He/H 

Thick-Target 

JS 

Po(MV) 

0 

2.7 - 3.3 

80 - 140 

5x10" 3 

2.4 - 3.0 

110 - 220 

0 

Thin-Target 
1.7 - 2.1 

120 - 220 

5x10-5 

1.5 - 1.9 

170 - 330 
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TABLE 5 


Average Particle Spectra at the Sun for the 1972, August 4 Flare 


Energy Range 


Electrons 


0.02 MeV < E < 0.14 MeV 
0.14 MeV < E < 1 MeV 


1 MeV < E 


Protons 


2.5 MeV < E < 200 MeV 


nN(E) (particles MeV~ 1 cnT 3 ) 

2. lxlO 44 E (MeV)" 2 * 5 
2.9X10 43 E (MeV) -3 * B 
r 1. 3X10 43 E (MeV )~ 2 
[ 2.3xl0* 2 exp [ -E(MeV)/5 ] 

■* 1.6x10*6 E (MeV )“ 3 
, 1. 9X10 43 (dP/dE) exp[-P(MV)/170] 
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TABLE 6 


Positron Emitter 
// 

1. n r< 

12 N 

13 N 

14 0 
is 0 
19 Ne 


Principal Positron Emitters 


Half Life 


1.5xl6 -6 sec 


sec 


20,5 min 


0.011 sec 


10 min 


71 sec 


2.06 min 


■Maximum 



Positron 

Production 

Threshold 

Energy (MeV) 

Mode 

(MeV) 

53 

♦ 0 . 

292.3 


p+ 4 He^n + ... 

185 


TT + ^JU + + V 


1 /s 

J-O? 

.1 g » 1 rv _ . 

p-K“ % o * 

41o4 


pP 4 N-P°C+. . . 

17.1 


p 4i2 C -> 10 c+... 

34.4 

0.92 

p-P 6 CK 11 C+. .. 

27.5 


pPf N-P 1 C+. . „ 

3.1 


p-P a C-> ll C+... 

17.9 

16.4 

pP 2 C-P 2 N+n 

19.6 

1.19 

pP 6 0-P 3 N+» . „ 

5.5 


pP 4 NP 3 N+. .. 

9.0 

1.8(99.4%) 

p-P 6 0-* 14 0+. . „ 

30.7 

4.1( 0.6%) 

pP 4 NP 4 0+n 

6.3 

1.74 

pP 6 0-P 5 0+. . . 

14.3 


a P 2 C-> l5 0+n 

2.8 

2.2 

Qf +l6 0 _X9 Ne+n 

3.75 


17.4 sec 
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FIGURE CAPTIONS 


1* Neutron production cross sections 

2. Partial neutron production modes in the thin- target model with 

power-law spectra 

3. Partial neutron production modes in the thin-target model with 

exponential spectra 

4* Total neutron production in the thin and thick-target models with 
power -law and exponential spectra 

5o Probabilities for neutron escape, decay, and capture in the solar 
atmosphere (solid lines), and photon yields per neutron (dashed 
lines) for no 3 He in the photosphere 

6, Probabilities for neutron escape, decay, and capture in the solar 
atmosphere (solid line), and photon yields per neutron (dashed 
lines) for 3 He/H - 5x10“ 5 

7* Neutron production energy spectra in the thin- target model with 
power-law spectra 

8* Neutron production energy spectra in the thin-target model with 
exponential spectra 

9* The average, f ,2* 2 MeV photon yield per neutron in the thin and 

thick-target model, with power-law and exponential spectra, and 
with 3 He/H = 0 and 3 He/H = 5x10" 5 in the photosphere „ 

10, 12 C* and 16 0* production cross sections 

llo 12 C* and l6 0* partial production modes in the thin-target model 
with power-law spectra 

12* 12 C* and 16 0* partial production modes in the thin-target model 

with exponential spectra 

13o Ratios of the 12 C* 4 * 43 yield to the total neutron yield for the thin 
and thick-target models, and power-law and exponential spectra. 

14. 7 Li production cross section in eta reactions 

15. Instantaneous bremsst rah lung emission from electrons with power- 

law spectra 

16 a Ratios of the 4.4 MeV line intensity. to the 2.2 MeV line intensity 

in the thin and thick- target models, with power-law and exponential 
spectra, and with 3 He/H = 0 in the photosphere. The shaded region 
is the data of Chupp et al„ (1975) 
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17 o Ratios of the 4.4 MeV line intensity to the 2.2 MeV line intensity 
in the thin and thick-target models, with power- law and exponential 
spectra, and with 3 He/H = 5x10" 6 in the photosphere. The shaded 
area is the data of Chupp et al, (1975).. 

18 o The 2.2 MeV line intensity at Earth in the thin and thick-target 
models for power-law and exponential spectra and 3 He/H =* 0. 

19. The 2. 2 MeV line intensity at Earth in the thin and thick-target 

models for power-law and exponential spectra and 3 He/H = 5xl0“ 5 

20. The observed hard x-ray and gamma-ray continuum emission from 

the 1972, August 4 flare. The shaded area is data from van Seek 
et al. (1973), and the data points are from Suri et al e (1974). 

Curves A, B and C are the bremsst rah lung spectra produced by the 
corresponding electron spectra shown in Figure 21. 

21. Instantaneous electron and proton numbers at the Sun* The electron 

spectra A, B and C (solid lines) produce the corresponding bremsstrah- 
lung spectra in Figure 20. The dashed lines are proton numbers 
deduced from line emissions in the thin-target model with 
3 He/H = 0. 

22. Energy loss per 4.4 MeV photon at Earth in the thin and thick-target 

models with power-law and exponential spectra* In both models 
the power-law spectra are of the form kE“ s for E>Ec and kEc~ s for 
E<Ec, 

23 * Radioactive positron emitter production cross sections* 

24* Partial positron production modes in the thin-target model with 
power-law spectra 

25* Ratios of the tt* meson yields to the total positron yields for the 
thin and thick-target models, and the power-law and exponential 
spectra 

26 „ Ratios of the total positron yield to the total neutron yield for 
the thin and thick-target models, and power law and exponential 
spectra 

27. Positron free -annihilation rate and positronium formation rate in 
a hydrogen plasma 

28* Instantaneous pair production rates by relativistic electrons 

29. Ratios of the gamma-ray yield from TT° decay to the total neutron 

yield for the thin and thick-target models and power- law and expon- 
ential spectra 

30 o Ratios at Earth of the time -integrated fluxes of neutrons to 2*2 MeV 
photons, in the thin and thick- target models with power- law and 
exponential spectra for isotopic neutron release at the Sun and 
3 He/H = 5xlCT 5 
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